Intermittency route to chaos for the nuclear billiard - a quantitative study 
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We extended a previous qualitative study of the intermittent behaviour of a chaotical nucleonic 
system, by adding a few quantitative analyses: of the configuration and kinetic energy spaces, power 
spectra, Shannon entropies, and Lyapunov exponents. The system is regarded as a classical "nuclear 
billiard" with an oscillating surface of a 2D Woods-Saxon potential well. For the monopole and 
dipole vibrational modes we bring new arguments in favour of the idea that the degree of chaoticity 
increases when shifting the oscillation frequency from the adiabatic to the resonance stage of the 
interaction. The order-chaos-order-chaos sequence is also thoroughly investigated and we find that, 
for the monopole deformation case, an intermittency pattern is again found. Moreover, coupling 
between one-nucleon and collective degrees of freedom is proved to be essential in obtaining chaotic 
states. 

PACS numbers: 24.60.Lz, 05.45.-a, 05.45.Pq, 21.10.Re 
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I. INTRODUCTION 

We begin by briefly reminding that a conjugated con- 
tinuous effort has been made to relate the emergence of 
the collective energy dissipation through one and two- 
body nuclear processes with the chaotical behaviour of 
nuclear systems [H-[34j|. 

A few options in choosing the collective oscillation 
frequencies have come into focus in the past years, in 
connection with the onset of chaoticity for "nuclear bil- 
liards" . First of all, the issue of dissipation into thermal 
motion of the adiabatic collective vibrational energy of 
the potential well was treated for several multipolarities 
by Burgio, Baldo et al. [J-Q- On the other hand, when 
trying to associate different vibration frequencies to var- 
ious nuclear processes, the path to chaos was found to be 
changed with the order of multipole (29l-[34j . 

This paper was intended to bring a quantitative ar- 
gumentation, based on a systematic study of the con- 
figuration and kinetic energy spaces, power spectra, in- 
formational entropies, and largest Lyapunov exponents. 
The study was done in completion of a few qualitative 
types of analysis previously presented [34J: sensitive de- 
pendence on the initial conditions, single-particle phase 
space maps, fractal dimensions of Poincare maps, and 
autocorrelation functions. 

In short, we remind that, by studying the nucleonic 
dynamics in a Woods-Saxon potential, one can find an 
increase of the chaotical degree of the system behaviour 
as raising the frequency of 2D wall oscillation. 



The main result of [341 ] was reported in relation with 
an intermittent route to chaos for the monopole vibra- 
tions close to the resonance phase of a nuclear interaction. 
Still, we mention that the purpose of these two coupled 
etudes was only to emphasize the detection of such inter- 
mission for the "nuclear billiards" and not to establish 
its type according to [3a], nor to compare it with other 
intermittency patterns from known experimental results 



II. TOY MODEL 

We continue the study on a classical dynamical system 
proposed by Burgio, Baldo et al. [IH1], system composed 
of a number of A nucleons with no charge, spin, or in- 
ternal structure. A two-dimensional deep Woods-Saxon 
potential well, regarded as " nuclear billiard" , is periodi- 
cally hit by the nucleons. The Bohr Hamiltonian in polar 
coordinates is a sum of two components: kinetic (Ekin.) 
and potential (E po t.), the kinetic one decoupling into ra- 
dial (E r ), centrifugal (El), and collective terms (E co u): 



Ekin. — E r + El + E co iL — 



^— ' \ 2m 2mr 2 - 



2M' 
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Epot =Yv(r h R{e j )) + ^^. (2) 
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The phase space is defined by particle and collective 
momenta and their conjugate coordinates: (r,p r ), (9,pe) 
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and (a, p a ). The collective coordinate a oscillates with 
f2 frequency, the Inglis mass M is equal to itiARq, and 
the nucleon mass: m = 938 MeV. 

For the time being we are not interested in studying 
the nucleon dynamics beyond the Woods-Saxon barrier: 



F(r i)J R(%)) = 



V 



1 + exp 



ri-R(8j,a) 



(3) 



and therefore we choose a deep well: Vo = —1500 MeV 
and accordingly, a low value for the diffusivity coefficient: 
a = 0.01 fm. 

When considering the two-dimensional case, the fron- 
tier of the collective motion is described as a function of 
the collective variable a and of the Legendre polynomials 
Pl (cos Bj) 018: 



Rj = R (6 3 , a) = Rq [1 + aP L (cos Oj)] 



(4) 



The oscillation degree of the potential well L is con- 
sidered for the monopole (0), dipole (1) and quadrupole 
case (2). 

If the surface has a stationary behaviour, or whenever 
one takes into account the uncoupled Hamilton equations 
(UCE) for the particle: 



Pr, 



Pn 



dV_ 
r j ' 



dV dRi 
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n or-j 



dRj d8j 



and collective degrees of freedom (d.o.f.): 
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dV 
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da 



(5) 



(G) 



Rq has a fix value, chosen for consistency with previous 
papers [H-d, as 6 fm. 

A Runge-Kutta type algorithm (order 2-3) with an op- 
timized step size was used for solving the system of dif- 
ferential equations, while keeping the absolute errors for 
the phase space variables under 10 -6 and conserving the 
total energy with relative error: AE/E w 10~ 8 (Fig. [l}. 

We imposed the equilibrium condition between the 
pressure exerted by the particles and the mechanical pres- 
sure of the wall [l|-y, [34| and thus obtained for the ini- 
tial equilibrium value of the collective variable, perturbed 
with a small value: 



a 



-1 + yj\ + 8T/mRlfl 2 



0.15, 



(7) 



where T = 36 MeV [l-|3j, l3J is the two-dimensional 
kinetic energy and also the temperature of the nuclear 
system, when considering the natural system of units 
(h = c = kg = !)■ 
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FIG. 1. Energy conservation for the adiabatic and resonance 
phase of the interaction (L = 0). 



III. THE QUANTITATIVE ANALYSIS OF THE 
ROUTE TO CHAOS 



We carry on the study begun in [34j by gradually 
changing the degree of vibration of the potential wall 
from a slow motion (adiabatic state) to a rapid one (the 
so-called resonance state of the interaction). In the first 
case, the collective motion is described by a radian fre- 
quency smaller than 0.05 c/fm. The latter is dominated 
by frequencies close to the one-nucleon collisional fre- 
quency: 



Ro 



Upart — — ■ y — ~ 0.145 c/fm. 



(8) 



The physical motivation for studying the one-nucleon 
chaotical dynamics in the " nuclear billiard" , ranging the 
frequencies from the adiabatic to the resonance regime of 
a nuclear interaction, is explained in some detail in [34j ] . 

We briefly remind that the process of nuclear multi- 
fragmentation can be viewed as a resonance process and 
that for smaller excitation energies, nuclear evaporation 
or breakup of a projectile nucleus occurs when the energy 
is shared between the collective and one-nucleon degrees 
of freedom. 

By using a few types of analyses: sensitive depen- 
dence on the initial conditions, single-particle phase 
space maps, fractal dimensions of Poincare maps and au- 
tocorrelation functions, we emphasized that an intermit- 
tent route to chaos is observed in the monopole case when 
increasing the vibrational frequency to Q = 0.1 c/fm [34| . 
In the resonance phase of the interaction the onset of 
chaotical behaviour was found to be earlier than at any 
other adiabatic oscillations of the Woods-Saxon potential 
well. 

We present here other methods [46j promoting the idea 
that the degree of chaoticity increases when moving from 
the adiabatic to the resonance regime: analyses of the 
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configuration and kinetic energy spaces, power spectra, 
generalized informational entropies and Lyapunov expo- 
nents. Furthermore, we try to identify possible pathways 
to chaos, including the intermittent one, previously put 
in evidence [341 ]. 

A. Configuration and kinetic energy spaces 

In order to establish if a specific physical system 
presents a chaotic dynamics and to identify possible 
routes to chaos we analyzed the behaviour of a small 
bunch of trajectories in the configuration and in the ki- 
netic energy space, respectively. For example, we took 
five trajectories separated by an e = Ar = 0.01 fm aper- 
ture, while keeping the rest of the initial phase space 
variables constant and let the system evolve over a given 
time (At = 1,600 fm/c). 

For the transient stages from adiabatic to resonance, 
the temporal evolution of an initially confined trajectory 
bundle was studied for the monopole and dipole oscilla- 
tion modes of the potential wall and also for the limit 
situation, in which the individual and collective degrees 
of freedom remain uncoupled (Figs. and [3]). 

The configuration space revealed a high degree of sym- 
metry in (r, 9) plane in both cases, (4 + 2) uncoupled 
nonlinear differential equations and monopole (left and 
middle panels). Also, that the central zone remained 
uncovered, reflecting the conservation of the nucleon an- 
gular momentum. 

Another important conclusion was issued from the defi- 
nition of the stability concept of a dynamical system. For 
the aforementioned cases the dispersed trajectory pack 
periodically regroups on the frontier that delimits the for- 
bidden zone of the phase space. This type of behaviour 
corresponds with the definition of stability given by Pois- 
son (for e.g., in ?47j). We will herewith remind that the 
Poisson stability defines as steady the movement of a 
particle system of which configuration comes close, from 
time to time, to the initial position. 

At a first glance, on the simple UCE case one can dis- 
tinguish two extremities of the radius of the particle pe- 
riodic motion in the 2D potential well: r m i n (UCE) w 
1.42 fm, and rMax(UCE) « 5.96 fm. These values cor- 
respond to the roots of the boundary equation: 



for a given one-nucleon energy E, when the radial com- 
ponent of velocity vanishes (r= 0). 

The analysis of the particle motion in the configura- 
tion space is similar to that applied to any system with 
bound unclosed trajectories. The nucleonic motion takes 
place within a circular crown (the so-called annulus) de- 
termined by the concentric circles of r m i n and ruax radii. 
The trajectory is symmetric about any turning point. 



For At = 1,600 fm/c we descry in the UCE case as 
much as 27 distinct apocenters (at r — rj\j ax ). These 
are correlated with a number of 13 complete and one 
incomplete revolutions about the center of the force field 
(i.e. 27 straight lines before r(t) changes its sense of 
variation). 

The particle completely sweeps over twice the 2D con- 
figuration space after 1,537 fm/c. However, we notice 
that the bound trajectories are open, which means that 
the orbits never pass twice through a given point (see Fig. 
[2]), which is in concordance with Bertrand's theorem. We 
briefly remind that for a bound orbit to be closed, the an- 
gle between two consecutive apocenters must be: 

Ad = 2ir-—, (10) 
n 2 

i.e. after n 2 revolutions about the center, the radius 
vector should sweep out a multiple ri\ of 2-7T radians. For 
the UCE case above considered we consequently obtain: 
AOuce ~ 47r/13 radians (Fig. [2]- left column). 

The kinetic energy points are displayed in right isosce- 
les triangular shaped patterns (Fig. [3]), whose hy- 
potenuses are described by Eqs. (11) and (12), for the 
two specific non-chaotical situations: UCE and the inter- 
mittent monopolar "window" emerged at Q = 0.1 c/fm: 

E LuCE = 18.05 - E r , (11) 

£Wic /fm = 18.09 -£ r . (12) 

Moreover, for the intermittency frequency of monopole 
oscillations, one can notice a second smaller segment with 
the same negative slope: 

E Ln=0 . lc/tia =5.87-E r . (13) 

For the uncoupled differential equations there are a 
couple of extreme values for the centrifugal kinetic en- 
ergy: E Lmin (UCE) = 1.02 MeV, and E Luax {UCE) = 
18.05 MeV, associated with tmo,x{UCE) and respec- 
tively, with r m i n (UCE) (left column of Fig. [3] and Eq. 

As for the monopolar intermittency, we can distinguish 
just five distinct values for the E L : 1.32 MeV, 2.13 MeV, 
3.01 MeV, 5.87 MeV, and 18.09 MeV (central plot of Fig. 
[3]), correlated with stationary radii: r\ ~ 5.26 fm, r 2 ~ 
4.14 fm, r 3 w 3.49 fm r 4 « 2.51 fm, and r 5 » 1.43 fm 
(Eq. (1), Fig. 2 of [33, and central plot of Fig. |2J). 
Thus, the nucleonic motion for the intermittent case is 
composed of alternated revolutions about the force field 
centre, forming a cyclic symmetrical structure, for e.g., 
ri, r 5 , r 2 , r 4 , r 3 , r 4 , r 2 , r 5 , r 4 , and so on (Fig. [2]). This 
behaviour can be easily verified through the sensitivity 
dependence on the initial conditions analysis, previously 
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FIG. 2. The variance in time (At — 1, 600 fm/c) of a packet of one-particle trajectories in the configuration space for adiabatic 
(upper row), intermediate (middle row), and resonance (lower row) collective vibrations for the uncoupled single and collective 
motion, monopole, and dipole cases. 



presented (third column of Fig. 2 - (34[). It should also 
be mentioned that, following this radius alternation, the 
nucleon covers in 2D configuration space « 2n radians 
after 6 full revolutions in almost 590 fm/c. 

Concluding, we highlight once more, that ordered, non- 
chaotical events, exhibit periodical symmetrical patterns 
in the configuration and kinetic energy spaces. This was 
shown to be a characteristic feature of the uncoupled non- 
linear Hamilton equations case and also, of the steady, 
intermittent behaviour arisen in the monopole case at 
0.1 c/fm vibrational radian frequency. 

A tendency to compactly fill the kinetic energy 
space when increasing the monopolar vibrations (from 
0.05 c/fm to 0.145 c/fm) was observed, except for the 
intermittency situation above described. For the L = 1 
oscillation mode of the potential well, it seems that at the 
same frequency (51 = 0.1 c/fm) a somewhat intermittent 
behaviour could also come out, but this was proved to be 
elusive, as verified when reverting to this issue with the 
help of informational entropies and Lyapunov exponents 



and analyzing the system on longer time periods. 

B. Power spectra 

In order to better distinguish between a multiple peri- 
odical behaviour that can also exhibit an erratic pattern 
and chaos we used the Fourier transform of the analyzed 
signals: 

x (u) = lim / e luJt ■ x (t) dt , (14) 

N 

x (lu) = lim V e iuJt " ■ x [t n ) . (15) 

71 = 

For a multiple periodical movement the power spec- 
trum: 
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FIG. 3. The variance in time (At — 1, 600 fm/c) of a packet of one-particle trajectories in the kinetic energy space for adiabatic 
(upper), intermediate (middle), and resonance (lower) collective vibrations for the uncoupled single and collective motion, 
monopole, and dipole cases. 



P{oj) = \x(uj)\ 



(16) 



will only contain a number of discrete lines: the fun- 
damental frequencies of the system and their associated 
sets of harmonics, while the chaotical behaviour is com- 
pletely aperiodical and is represented by a continuous or 
quasi-continuous broadband. 

The obtained results are presented in Figures |4][7] As 
a persistent feature of the physical system analyzed one 
should mention that for the monopole and dipole defor- 
mation degrees of the potential well (Figures [5] and O 
the chaotic behaviour increases in time, thus confirming 
previous results. 

The transition towards a chaotic regime was put again 
in evidence once passing from the adiabatic to the reso- 
nance stage of the interaction. The power spectra reveal, 
as expected, the intermittent feature of the transition in 
the monopolar case at ft = 0.1 c/fm. 

This can be detected for periods of time large enough 



(At > 1,600 fm/c) to positively identify chaotic pat- 
terns, by transition from a quasi-continuous spectrum of 
the one-nucleon radial coordinate (il a d = 0.02 c/fm) to a 
discrete periodical one, containing fundamental frequen- 
cies of the system and its harmonics (Q = 0.1 c/fm) and 
again to a continuous spectrum at the resonance vibra- 
tional frequency (Q r es = 0.145 c/fm) (Fig. 0- right pan- 
els). The order-chaos-order-chaos sequence can be also 
spotted out for the monopole oscillations in the power 
spectra of the collective degree of freedom (Fig. 0- sec- 
ond column). 

The temporal series of the radius variable show a sym- 
metrical sawtooth waveform for the uncoupled situa- 
tion at any chosen vibration frequency, and also for the 
monopole case at adiabatic collective oscillations. For 
the rest, in general an asymmetrical sawtooth form de- 
fines the series, but sometimes, more complicated pat- 
terns appear at higher multipole orders (Figures 1-4 - 

LU). 

The difference between two successive maxima in the 
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FIG. 4. The power spectra of the one-nucleon radial variable for different phases of the interaction (from adiabatic to resonance 
regime) calculated for the uncoupled single and collective degrees of freedom case. 



temporal series of the radius variable for the UCE case 
is: T 0uCE w 1,537/26 = 59.1 fm/c (left column of Fig. 
[3]) and can be also obtained from the sensitive depen- 
dence on the initial conditions analysis (see Fig. 1 - 
[Hj]). The fundamental frequency for the radial saw- 
tooth temporal series: ujq uce — 2it/Tq uce = 0.106 c/fm 
and its first three harmonics: io\ UCE = 0.212 c/fm, 
uJ2 UC e = 0.318 c/fm, and (^3 UCE — 0.424 c/fm, can be 
easily traced down in Figured) 

As for the "window" of intermittency at L = 0, we 
obtained: T 0int ~ 590/12 = 49.2 fm/c (Fig. 2 - Q 
and central plot of current Fig. [5]). This gives the cor- 
responding fundamental frequency: u>Q int — 0.128 c/fm 
and its associated harmonics: u>i int = 0.256 c/fm, u>2 int — 
0.384 c/fm, and oj 3mt = 0.512 c/fm (Fig. |SJ). 



mational entropy (or Shannon entropy), introduced as 
usually [Il,[|l|53j]: 

N(t) 

Sshannon (*) = -^Pk ' lnPfc, (17) 
k=l 

N (t) being the number of gradually occupied cells un- 
til the time t. 

This type of entropy is actually a number which quan- 
tifies the time rate of information production for a chaotic 
trajectory [13]. We consider in the first place the case of 
a particle that at every moment occupies a cell of the 
two-dimensional lattice phase space with a pk probabil- 
ity: 



C. Shannon entropies 

In order to further investigate route to chaos, we paid where: 
attention to the time evolution of the generalized infor- 



Pk — 1/Ntotal cells, (18) 



Ntotal cells =N r -N Pr -N ei (19) 
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FIG. 5. The power spectra of the one-nucleon radial variable for different phases of the interaction (from adiabatic to resonance 
regime) calculated for the monopole case. 



N r , N Pr , and Ng are the number of bins of the (r, p r , 0) 
lattice. For pg is a constant of motion for the monopole 
and the UCE cases, we use for comparisons only these 
three phase space variables. 

As an alternative measure for the above defined en- 
tropy we also used the cumulative filling percentage of 
the one-nucleon phase space: 



V(t)= M W ■100(%). (20) 

+ *total cells 

In the first place, for a given wall frequency of vibra- 
tion and for a certain multipolarity (here, for £l res = 
0.145 c/fm and L = 0), we studied the dependence of the 
Shannon entropy with the number of bins. A clear ten- 
dency for smoothing the entropy curve was found when 
decreasing the bin. A reduced number of cells (N^ = 2 3 ) 
is characterized by an entropy formed from a small num- 
ber of high-amplitude Heaviside functions. As the num- 
ber of bins increases (for e.g., here to 12 3 ), the entropy 
gets a more realistic representation, being composed of 



a superior number of low-amplitude step functions (Fig. 

ED. 

Moreover, the filling percentage rj of the one-nucleon 
phase space maps can drastically differ with the size of 
the bin. Thus, after the system evolved over 400 fm/c, 
a phase space with 8 bins is entirely covered, 64 bins 
can be filled in with 0.8594 probability, a 26.95 filling 
percentage for 512 cells can be found, and we counted 
only as much as 226 bins occupied out of a total of 1, 728 
(i.e. n = 13.08 %). 

At a first glance one can identify a series of entropy 
plateaus, which could be put in correspondence with sta- 
tionary or quasi-stationary thermodynamic values of the 
system if a large number of particles would be under 
study. Some of them will vanish when considering a large 
number of bins. However, those surviving for Nf, — > cxo 
could be associated with stationary nucleonic states in 
the chosen potential well in the limit of a large number 
of degrees of freedom. 

For a given 2D phase space lattice formed of Nb = 4 3 
bins we present in Figures [911121 a comparison between 
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FIG. 6. The power spectra of the one-nucleon radial variable for different phases of the interaction (from adiabatic to resonance 
regime) calculated for the dipole case. 



the informational entropies of the physical system in 
study, starting from the adiabatic stage of interaction 
and gradually increasing the vibrational wall frequency 
towards the dipole resonance value, £l res = 0.145 c/fm. 
The slopes for the resonance frequency case were found 
to be significantly higher than for the adiabatic one 
(^ad = 0.02 c/fm) for all multipolarities involved. 

Another comparison revealed significant differences be- 
tween the onset times of the quasi-constant Shannon en- 
tropy values for all cases taken into consideration. Thus, 
for four vibrational radian frequencies and for four cou- 
pling modes of the Hamilton equations we show the infor- 
mational entropy values after 800 fm/c (Table HJ and the 
associated phase space filling degrees (Table HI)). Also, in 
Table IIII1 are presented the periods of time after which 
the filling percentages r\ equal unity. 

We continue the analysis by further defining the Shan- 
non entropy for a group of w nearby orbits: 



Straject. pack (t) — In N w (i) , 



(21) 



TABLE I. The computed S S k„„„™ (t = 800 fm/c) of the 
(r <-> p r «-» 9) one-particle phase space maps at several multi- 
polarities and frequencies of wall vibration 



Oscillation frequency 


UCE 


L = 


L = 1 


L = 2 


Q.ad = 0.020 c/fm 


3.6889 


3.7136 


3.7842 


3.4340 


Q ad = 0.050 c/fm 


3.6889 


4.0775 


4.0775 


3.2958 


n =0.100 c/fm 


3.6889 


3.6889 


3.8501 


4.0073 


fires = 0.145 c/fm 


3.6889 


4.1589 


4.0431 


4.0254 



so that the number of occupied cells is: 

1 < N w (t) < w, 



(22) 



thus describing the spread of the trajectories at each 
moment of time t (Figs. ITSHTB)) . When reaching the 
maximum divergence, the entropy for five distinct phase 
space paths gets its highest value (Table ITVT) . 
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Uncoupled Eqs. 



L = 



L= 1 



I i U . A a 



0.2 0.4 

x 10 7 



3 



1 

o- 

0.2 0.4 

x10 6 



10 

3, 

Q_ 5 



10 

Q_ 5 



0.2 0.4 

x10 6 



i 1 1 L 



0.2 0.4 

(o (c/fm) 




L = 2 



0.2 0.4 

(o (c/fm) 



0.2 0.4 

(o (c/fm) 



x 10 



0.2 0.4 

x10 e 



il 



0.2 0.4 

x10 e 



0.2 0.4 




0.2 0.4 

(o (c/fm) 



FIG. 7. The power spectra of the collective coordinate for different phases of the interaction (from adiabatic to resonance 
regime) calculated for the uncoupled single and collective degrees of freedom, monopole, dipole and quadrupole cases. 




FIG. 8. The Shannon entropy (upper panel) and the filling 
degree of the one-particle phase space (lower panel) for all 
frequencies studied (fl = 0.02 — 0.145 c/fm), computed for 
the uncoupled single and collective degrees of freedom case. 
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FIG. 9. The Shannon entropy (upper panel) and the filling 
degree of the one-particle phase space (lower panel) for all 
frequencies studied (Q — 0.02 — 0.145 c/fm), computed for 
the uncoupled single and collective degrees of freedom case. 
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FIG. 10. The Shannon entropy (upper panel) and the filling 
degree of the one-particle phase space (lower panel) for all 
frequencies studied (Q = 0.02 — 0.145 c/fm), computed for 
the first degree of multipolarity (L = 0) . 
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FIG. 11. The Shannon entropy (upper panel) and the filling 
degree of the one-particle phase space (lower panel) for all 
frequencies studied (Q = 0.02 — 0.145 c/fm), computed for 
the dipole case (L = 1). 



TABLE II. The filling percentage rj of the (r <-> p r <-> 9) one- 
particle phase space maps at several multipolarities and fre- 
quencies of wall vibration 



Oscillation frequency 


UCE 


L = 


L = 1 


L = 2 


fi ad = 0.020 c/fm 


62.50 


64.06 


68.75 


48.44 


fi ad = 0.050 c/fm 


62.50 


92.19 


92.19 


42.19 


O =0.100 c/fm 


62.50 


62.50 


73.44 


85.94 


n res = 0.145 c/fm 


62.50 


100.00 


89.06 


87.50 




t (fm/c) 



FIG. 12. The Shannon entropy (upper panel) and the filling 
degree of the one-particle phase space (lower panel) for all 
frequencies studied (fi = 0.02 — 0.145 c/fm), computed for 
the quadrupole oscillations (L = 2) . 



TABLE III. The time (in fm/c) at which the informational 
entropies of the (r <-> p r <-> 9) one-particle phase space maps 
at several multipolarities and frequencies of wall vibration 
have the maximum value (i.e. r\ = 100 %) 



Oscillation frequency 


UCE 


L = 


L = 1 


L = 2 


fi ad = 0.020 c/fm 


> 10 s 


6,023 


6,359 


5,356 


Q ad = 0.050 c/fm 


> 10 s 


1,618 


4,223 


> 10 5 


n =0.100 c/fm 


> 10 8 


11,442 


3,241 


2,758 


fires = 0.145 c/fm 


> 10 5 


729 


1,887 


10,571 





t (fm/c) 



FIG. 13. The Shannon entropy of a bunch of five one-particle 
close trajectories for frequencies between 0.02 and 0.145 c/fm 
(uncoupled single and collective degrees of freedom case). 
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FIG. 14. The Shannon entropy of a bunch of five one-particle 
close trajectories for frequencies between 0.02 and 0.145 c/fm 
(L = 0). 



FIG. 16. The Shannon entropy of a bunch of five one-particle 
close trajectories for frequencies between 0.02 and 0.145 c/fm 
(L = 2). 
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FIG. 15. The Shannon entropy of a bunch of five one-particle 
close trajectories for frequencies between 0.02 and 0.145 c/fm 
(L = 1). 



TABLE IV. The time (in fm/c) at which the one-particle 
Shannon entropies of a pack of w = 5 close orbits begin having 
the maximum value (i.e. Straject. pack = 1.60944) for various 
coupling degrees between the one-nucleon and the collective 
d.o.f. and for the standard chosen wall frequencies 



Oscillation frequency 


UCE 


L = 


L = 1 


L = 2 


Q ad = 0.020 c/fm 


> 10 4 


1,095 


555 


688 


fi ad = 0.050 c/fm 


> 10 4 


855 


476 


122 


fi =0.100 c/fm 


> 10 4 


4,133 


333 


395 


fi res = 0.145 c/fm 


> 10 4 


279 


327 


469 



We begin the analysis with the UCE case. The single 
and collective uncoupled d.o.f. give birth to a quasi- 
laminar behaviour with a weak development of chaotic 
states. The one-particle informational entropy shows an 
identical evolution, no matter the frequency chosen. The 
orbit covers, after 800 fm/c, only 62.50 % of the entire 
lattice (Table HIl and Figure[9]) and does not reach 100 %, 
even after At = 100,000 fm/c (Table |m|- Also, the 
phase space is not covered up by all five trajectories for 
the whole range of 10, 000 fm/c considered, when analyz- 
ing Straject. pack 

(Fig. [I3]and Table [IV}. 

For the dipole oscillations mode, at Q, a d = 0.05 c/fm, 
it appears that, after only 800 fm/c, the entropy closes 
in upon its maximum value: S M ax = In N total ceils = 
4.1589 (Fig. [TT1 and Table Q|. However, on long periods 
of time, the real tendency is towards filling up the nu- 
cleonic phase space as rapid as the vibrational frequency 
is increased (Table HIT)) . The exact pattern is repeated 
when studying the Shannon entropy for closeby nucle- 
onic trajectories (Fig. [TBI and Table HV*]) . 

We found quite the same feature for the monopole 
case, with exception for the intermittent "window" at 
= 0.1 c/fm (Tables HI H] and Fig. [TDJ). The occupy- 
ing rate is so small in the intermittent zone, that just at 
11,442 fm/c, the particle would have covered the whole 
phase space (see Table lHTj) . A similar conclusion can be 
drawn from Table IfVI and Fig. [H] (with a double tem- 
poral scale scanned for the intermittent frequency). The 
trajectory pack informational entropy reaches its highest 
value after the longest one-particle evolution time of all: 
4, 133 fm/c. 

The quadrupole oscillation also reveals an apparent in- 
termittent pattern, this time at £l a d — 0.05 c/fm. We call 
it intermittent because after 800 fm/c the nucleon fills in 
only 42.19 % of the total number of bins (Figure [T2l and 
Table [TT]) . and a longer time than 100,000 fm/c is re- 
quired to get to r\ = 100 % (Table lilT|) . However, this 
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behaviour can be a misleading one, the Shannon entropy 
for a trajectory bunch showing exactly the opposite (see 
Fig. [Hand Table Hv). after 122 fm/c the orbits being 
completely dispersed. 



TABLE V. The largest Lyapunov exponents (in c/fm) com- 
puted as slopes of the increasing branches of the ordinary 
Lyapunov exponents for the degrees of multipole and vibra- 
tional frequencies used 



D. Lyapunov exponents 

We furthermore presented another quantitative analy- 
sis: the temporal evolution of the Lyapunov exponents, 
A (t). As previously shown, initial adjacent points in the 
phase space Axo(t — 0), can generate in time separated 
trajectories Ax(t). When studying the evolution of a 
single phase space parameter, the one-dimensional Lya- 
punov exponent takes the form: 



Oscillation frequency 







A(t) 



lim In 

\Ax o \-*0 



Ax (t) 



Ax 



(23) 



The generalization for obtaining the multi-dimensional 
Lyapunov exponent is then straightforward: 



X(t) = In 



EfcLi [xk (*) - Xko]' 
rloi 



(24) 



where the sum is taken over all to = 4 squared dif- 
ferences between final Xk (t) and initial Xko one-nucleon 
phase space variables. Integration times of the order of 
10 3 fm/c exclude errors when computing the Lyapunov 
exponents. 

In short, we here remind that the trajectories can be 
classified as function of the Lyapunov exponents. Thus, 
one can distinguish periodical behaviours, for A = 0, dis- 
sipative movements with a hxed point or a basin of at- 
traction (A < 0), and aperiodical chaotic states (A > 0), 
when the iterative discrete evolution of the solution series 
(Eqs. 5 and 6) leads to a chaotic pattern. 

Another way of measuring the system sensitivity to 
initial conditions is to compute the largest Lyapunov ex- 
ponent (LLE). Usually a couple of methods can be em- 
ployed, one based on the time dependence of the multi- 
dimensional Lyapunov exponent, the other on Wolf's 
standard method that uses a Gram-Schmidt Reorthonor- 
malization of the tangent vectors [Hj]. In the latter, the 
LLE is obtained by taking the asymptotic value of the 
multi-dimensional Lyapunov exponent: 



fi ad = 0.020 c/fm 
n ad = 0.050 c/fm 
n =0.100 c/fm 
ft res = 0.145 c/fm 



0.003939 
0.004432 
0.000761 
0.008739 



0.008689 
0.009829 
0.015454 
0.016662 



0.004306 
0.023203 
0.014402 
0.010086 



a certain time t c (Fig. [T7|) . The straight lines represent 
fits whose slopes match the LLE (Table |Vj) . They are in 
inverse proportion with the onset times of chaoticity (r = 
I/Xmux), being a measure of the trajectory decoupling at 
a microscopic level. 

When the single-particle and collective d.o.f. remain 
uncoupled, the Lyapunov exponents basically oscillate 
between two quasi-stationary regimes. This happens for 
all vibrational frequencies involved, reflecting a periodi- 
cal regrouping of orbits in two basins of attraction. The 
phase space not being covered, even after a hundred of 
thousand of fm/c, computing the LLE becomes futile for 
this case. 

One can remark for dipole oscillations (Fig. [17]- middle 
column) a faster evolution towards reaching saturation 
states of the 4-dimensional Lypaunov exponents, once 
passing from the adiabatic (r a( j = 115 fm/c) to the reso- 
nance phase of the interaction (r, es = 60 fm/c). 

In the monopolar case the intermittency can be eas- 
ily traced at 0.1 c/fm vibrational frequency (Fig. [T7] 
- left panels). During the intermittent stage, indepen- 
dent nearby orbits microscopically diverge with the slow- 
est rate of all: r = 1,314 fm/c (Table EJ. In order 
to catch the heaving in sight of the stationary plateau 
at w 4,917 fm/c, the temporal scale was scanned over 
6,400 fm/c. 

The study of the quadrupole collective oscillation case 
confirms the results obtained with all previous analyses. 
Namely, the neighbouring trajectories deviate one from 
each other after just 43 fm/c at an adiabatic frequency: 
0.05 c/fm. Also, when increasing f2, the LLE evolution 
pattern exactly matches that found with informational 
entropy measured for a group of orbits (Tables IIVI and 
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Xhiax = lim — — . (25) 

t—>oo t 

Still, this method has the disadvantage that the inte- 
gration times have to be at least an order of magnitude 
larger than those here considered. Other methods are 
slightly less efficacious, being more CPU time-consuming 
when simulating strong chaotic systems [55j | . 

We consequently used the first method and noticed the 
saturation behaviour, i.e. the arising of a plateau after 



IV. CONCLUSIONS 

We investigated the chaotic nucleonic behaviour in 
a two-dimensional deep Woods-Saxon potential well for 
specific phases of the nuclear interaction. By compar- 
ing the order-to-chaos transition for these cases of in- 
terest, from adiabatic to resonance regime, it was shown 
that the couplings between the one-particle dynamics and 
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FIG. 17. The temporal evolution of the Lyapunov exponents of the single-particle d.o.f. for multiple radian frequencies fi and 
multipolarities L of the nuclear surface. The LLE are given by fits (thickened solid lines) of the increasing portions of the 
4-dimensional Lyapunov exponents. 



high multipole vibrational modes significantly decrease 
the onset of the chaotic nucleonic motion towards realis- 
tic nuclear interaction time scale. 

The quantitative study enfolded a plethora of analyses, 
pointing out that the paths to chaos for the " nuclear bil- 
liard" are dissimilar for the studied multipolarities. For 
the first two multipole degrees we noticed a more rapid 
emergence of chaotic states as moving on towards higher 
radian frequencies of oscillation. When analyzing the sys- 
tem with quadrupolc collective deformations of the po- 
tential well, an order-strong chaos-weak chaos-order se- 
quence is revealed. Still, as emphasized in the " Shannon 
entropies" subsection, the quadrupole case represents an 
intricate one, and further analysis would be required be- 
fore concluding it. 

Every type of quantitative analysis strengthened previ- 
ous results regarding the monopolar intcrmittency route 
to chaos for the "nuclear billiard". The collective os- 



cillation frequency for the intermittent behaviour was 
located prior to the resonance state of interaction (at 
n = 0.1 c/fm). 

Further studies along the above issues are currently in 
progress. The used formalism can be improved by adding 
spin and charge to the nucleons. A semi-quantal treat- 
ment of this problem, including Pauli blocking effect, is 
hoped to shed more light on the discussed issue in the 
near future. 
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